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Ref.: 55600 SE 

Applicant : NordNav Technologies AB 
Spread Spectrum Signal Processing 

THE BACKGROUND OF THE INVENTION AND PRIOR ART 

The present invention relates generally to processing of spread 
spectrum signals by means of vector based algorithms. More 
particularly the invention relates to a method of receiving spread 
spectrum signals according to the preamble of clajm 1 and a 
signal receiver according to the preamble of claim 28. The 
invention also relates to a computer program according to claim 
26 and a computer readable medium according to claim 27. 

Spread spectrum technology is becoming increasingly important 
in many areas, of which cellular communication systems and 
global navigation satellite systems (GNSS) represent two 
important examples. Moreover, the large variety of transmission 
standards creates a demand for hybrid- or general-purpose 
receivers. For instance, both the cellular standards cdma2000 
and WCDMA involve transmission of spread spectrum signals. 
However, due to differences in the air interfaces one and the 
same terminal cannot be used in the two systems. Instead, a 
dedicated terminal must be employed in each system. Alter- 
natively, a dual mode terminal must be designed, which includes 
two separate transceiver chains. Correspondingly, a navigation 
receiver adapted for one GNSS, say the Global Positioning 
System (GPS; U.S. Government), is not able to receive signals 
from a satellite that belongs to a different GNSS, such as the 
Galileo system (the European programme for global navigation 
services) or the Global Orbiting Navigation Satellite System 
(GLONASS; Russian Federation Ministry of Defense). 



In order to make this possible, a multi-mode receiver must be 
designed. However, including multiple receiver chains in a 
single device is not only expensive, it also renders the unit bulky 
and heavy, particularly if more than two signal formats are to be 
processed. Therefore, a programmable software receiver 
solution is desired whose signal processing principles may be 
altered according to which signals that shall be received and 
processed. 

Various software solutions are already known for processing 
GNSS signals. The patent document WO02/50661 is one 
example in which a GPS receiver tracking system for guiding 
missiles is described. The receiving station here includes one or 
more processing stages for receiving and processing GPS data. 
A control signal controls whether a quantity of processing stages 
should be increased or decreased in response to the result of a 
correlation operation. The tracking and recovery of timing in- 
formation is entirely performed in software. However, the pro- 
cessing involves real-time calculation of relatively complex 
Fourier transforms, and therefore requires considerable pro- 
cessing resources. Thus, a receiver of this type is not parti- 
cularly suitable for small sized units, such cellular phones or 
portable GNSS-receivers. 

The article "TUPM 12.4: Software Solution of GPS Baseband 
Processing", ULSI Laboratory, Mitsubishi Electric Corporation, 
IEEE, 1998 by Asa! T., et al. outlines the principles for a soft- 
ware GPS receiver using an embedded microprocessor, which 
enables signals from up to eight satellites to be demodulated at 
around 40 MIPS (millions of / or Mega Instructions Per Second). 
A 2 MHz GPS baseband signal is here fed to the software radio 
receiver, which performs spectrum despreading and mixing; 
synchronous demodulation, satellite acquisition and tracking; as 
well as phase compensation. An algorithm is proposed through 
which the intermediate frequency (IF) signal is 1-bit sampled 
and where downconversion is performed before despreading. 
Resulting data in the form of 32 samples are continuously 
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processed in parallel. A DRAM (Dynamic Random Access 
Memory) stores a particular C/A (Coarse Acquisition) code table 
for each satellite from which signals are received. All other data 
and instructions are stored in a cache memory in order to 
5 accomplish a high-speed processing. Although this strategy is 
declared to result in a reasonably efficient implementation, the 
document lacks a specific description as to how the dispreading 
and the correlation should be effected in order to, in fact, obtain 
til is effect. 

10 Akos. D. et al., "Tuning In to GPS ~ Real-Time Software Radio 
Architectures for GPS Receivers", GPS World, July 2001 
describes a receiver architecture through which IF signal 
samples are fed directly from a radio front-end to a program- 
mable processor for continued processing. The article mentions 

15 the possibility of using single instruction multiple data (SIMD) 
instructions to process multiple data samples in parallel. How- 
ever, portability- and flexibility problems are recognized, and 
again, there is no teaching as to how an efficient parallel pro- 
cessing could actually be accomplished. 

20 Dovis, F. et al., "Design and Test-Bed Implementation of a 
Reconfigurable Receiver for Navigation Applications", Electro- 
nics Department, Politecnico di Torino, Navigation Signal Ana- 
lysis and Simulation Group, Spring of 2002 relates to the design 
of a reconfigurable GNSS receiver which is capable of fusing 

25 data from two or more different GNSS:s. The document sketches 
an architecture which, in addition to a radio front-end, includes a 
Field Programmable Gate Array (FPGA) and a Digital Signal 
Processor (DSP). The authors address various computational 
load issues. However, they do not present an algorithm which 

30 fulfills the identified requirements. 

Hence, the prior art includes a number of examples of software- 
based GNSS-receivers. Nevertheless, there is yet no distinct 
teaching of a highly efficient solution which is suitable for imple- 
mentation in software and has capabilities that are at least in par 
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with those of today's ASIC-based solutions for receiving and pro- 
cessing spread spectrum signals In real-time (ASIC = Application 
Specific Integrated Circuit). 

SUMMARY OF THE INVENTION 

5 The object of the present invention is therefore to provide a 
solution for receiving and processing spread spectrum signals, 
which solves the problems above and thus to present a truly 
software adapted strategy that may be implemented in a 
physically small and power efficient device. 

10 According to one aspect of the invention the object is achieved 
by a method of receiving spread spectrum signals as initially 
described, which is characterized by comprising a preparation 
for the correlation step, wherein the preparation takes place 
before the receiving the continuous signal. This preparation 

15 involves pre-generating a multitude of code vectors, which each 
represents a particular code sequence of the at least one signal 
source specific code sequence. Moreover, according to the 
invention, the correlation step involves multiplying at least each 
vector in a sub-group of the code vectors with at least one 

20 vector, which is derived from the data word. 

This strategy is advantageous because the proposed vector 
approach makes it possible for a digital processor (e.g. a micro- 
processor) to process multiple signal samples in parallel during 
each clock cycle and thus make very efficient use of the 
25 processor. Furthermore, by pre-generating the code vectors 
valuable processing capacity is saved, which In turn renders it 
possible to deal with a sufficient amount of input data per time 
unit in order to, for example track navigation signals that are 
transmitted by the satellites in a navigation satellite system. 

30 According to a preferred embodiment of this aspect of the 
invention, each code vector represents a particular signal 
source specific code sequence, which is sampled at the basic 
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sampling rate and quantised witli tiie quantising process that is 
used to produce the level-discrete sample values. Preferably, 
the code vectors are also further adapted to the format of the 
incoming data signal in a way which renders the following 
5 processing efficient. 

According to another preferred embodiment of this aspect of the 
invention, the at least one signal source specific code sequence 
represents so-called pseudo random noise. This type of code 
sequences is beneficial in environments where two or more 

10 different signals are modulated onto the same carrier frequency. 
Namely, the different pseudo random noise signals are ortho- 
gonal (or at least almost orthogonal) to each other. Conse- 
quently, a first informatioh signal spread by hieans of a first 
pseudo random noise signal causes a minimal deterioration of a 

15 second information signal which is spread by means of a second 
pseudo random noise signal, and vice versa. 

According to another preferred embodiment of this aspect of the 
invention, the receiving step involves down conversion of an 
incoming high-frequency signal to an intermediate frequency 

20 signal. The high-frequency signal is presumed to have a spec- 
trum which is symmetric around a first frequency and the inter- 
mediate frequency signal is presumed to have a spectrum which 
is symmetric around a second frequency, which is considerably 
lower than the first frequency. Hence, the receiving step 

25 transforms the information carrying signal to a baseband which 
is technically more easy to handle than the band at which the 
original incoming signal is located. 

According to another preferred embodiment of this aspect of the 
invention, the method includes the following steps. First, a 
30 maximum frequency variation of the second frequency due to 
Doppler-effects is determined. Then, a Doppler frequency interval 
around the second frequency is defined. The Doppler frequency 
interval has a lowest frequency limit equal to the difference 
between the second frequency and the maximum frequency 
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variation, and a highest frequency limit equal to the sum of the 
second frequency and the maximum frequency variation. After that, 
the Doppler frequency interval is divided into an integer number of 
equidistant frequency steps, and subsequently a frequency 
5 candidate vector is defined for each frequency step. This way of 
dividing the spectrum is advantageous because it provides a 
flexible modeling of the intermediate carrier frequency and its 
variations. 

According to another preferred embodiment of this aspect of the 
10 invention, an integer number of initial phase positions is 
determined for the i-ntormodiate frequency signa l. The integer 
number thus represents a phase resolution, such that a 
reTatively~hig¥ humbe 

phase resolution, whereas a relatively low number corresponds 
15 to a comparatively low phase resolution. Furthermore, a carrier 
frequency-phase candidate vector is defined for each combi- 
nation of carrier frequency candidate vector and initial phase 
position. The different frequency-phase candidate vectors are 
beneficial, since they allow the determination of a very accurate 
20 estimate of frequency- and phase characteristics of a received 
signal, which in turn vouches for a demodulated signal having a 
high quality. 

According to yet another preferred embodiment of this aspect of 
the invention, the number of elements in each carrier frequency- 
25 phase candidate vector is determined. Subsequently, the carrier 
frequency-phase candidate vectors are stored according to a 

I data format, which is adapted to performing correlation mWpka 
operations between the at le a s t u iil vllIui doiivud Fp grn the data 
word and a segment of a carrier frequency-phase candidate 
30 vector. This adaptation is advantageous because thereby 
several signal samples may be processed in parallel by means 
of comparatively simple operators. Preferably, the adapting of 
the data format involves adding at least one element to each 
segment of the carrier frequency-phase candidate vector, such 
35 that the segment attains a number of elements which is equal to 
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the number of elements in each of the at least one vector that is 
derived from the data word. Hence, it is namely possible to 
process a segment of the carrier frequency-phase candidate 
vector together with one of the at least one vector by either a 
5 so-called Single Instruction Multiple Data (S[MD)-operation or a 
logical XOR-operation. 

According to still another preferred embodiment of this aspect of 
the invention, the method also includes the following steps. 
First, a maximum variation of trie code rate due to Doppler- 
10 effects is determined. Then, a Doppler rate interval is defined 
around a center code rate. The Doppler frequency interval has a 
lowest code rate limit equal to the difference between the center 
coae"rlite~aha~tlie "ma^^^ 

frequency limit equal to the sum of the center code rate and the 
15 maximum code rate variation. In analogy with the intermediate 
frequency spectrum, the Doppler rate interval is divided into an 
integer number of equidistant code rate steps, and a code rate 
candidate Is defined for each code rate step. In further analogy 
with the intermediate frequency, an integer number of possible 
20 initial code phase positions is determined for each code rate 
candidate. The integer number thus represents a code phase 
resolution, such that a relatively high number corresponds to a 
comparatively high code phase resolution, whereas a relatively 
low number corresponds to a comparatively low code phase 
25 resolution. Typically, the spread spectrum code modulation is 
less sensitive to distortion caused by Doppler effects than the 
frequency modulation. Therefore, the code rate steps may be 
relatively large. However, a certain degree of Doppler shift esti- 
mation is desirable to attain a demodulated signal of high quality. 

30 According to yet another preferred embodiment of this aspect of 
the invention, a set of code vectors is generated for each signal 
source specific code sequence by sampling each code rate- 
phase candidate vector with the basic sampling rate, and thus 
produce a corresponding code vector. Again, this simplifies a 

35 parallel processing of several signal samples by means of 
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comparatively simple operators. 

According to still another preferred embodiment of this aspect of 
the invention, a modified code vector is generated on basis of 
each code vector by: copying a particular number of elements 
5 from the end of an original code vector to the beginning of the 
modified code vector, and copying the particular number of 
elements from the beginning of the original code vector to the 
end of the code vector. This extension of the code vector is 
highly advantageous because it allows extraction of tracking 
10 parameters through early-, prompt- and late techniques by 
performing a simple translation of local copies of the signal 
source specific code sequence. 

According to still another preferred embodiment of this aspect of 
the invention, a set of modified code vectors is stored for each 

15 signal source specific code sequence. Each modified code 
vector here contains a number of elements, which represent a 
sampled version of at least one full code sequence. As 
mentioned above, the modified code vector also includes a 
repetition of the beginning and the end of the sequence. A 

20 particular modified code vector is defined for each combination 
of code rate candidate and code phase position. These pre- 
generated modified code vectors save considerable processing 
capacity in the processor, since they need not be calculated in 
real time, but may simply be acquired from a memory. 

25 According to yet another preferred embodiment of this aspect of 
the invention, the data format of the modified code vectors is 
adapted with respect to the data format of the at least one 
vector that has been derived from the data word, such that a 
modified code vector and one of the at least one vector may be 

30 processed jointly by a SIMD-operation or an XOR-operation. 
Naturally, this is advantageous from a processing efficiency 
point-of-view. 

According to yet another preferred embodiment of this aspect of 
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the invention, the method involves an Initial acquisition phase 
and a subsequent tracking phase. The acquisition phase esta- 
blishes a set of preliminary parameters which are required for 
initiating a decoding of signals that are received during the 
6 tracking phafee. The parameters may include: a modified code 
vector, a carrier frequency candidate vector an initial phase 
position, a code phase position and code index, which denotes a 
starting sample value for the modified code vector. A successful 
acquisition phase thus results in that at least one signal source 
10 specific code sequence is identified and that this signal is 
possible to track by the receiver (i.e. may be continued to be 
received). 

According to stiri an^^^^ 

the Invention, the tracking in turn, involves the following steps. 

15 First, based the tracking characteristics, a prompt pointer is 
calculated for each modified code vector. The prompt pointer 
indicates a start position for the code sequence. The initial 
prompt pointer value is set equal to the code index. Then, at 
least one pair of early- and late pointers is assigned around 

20 each prompt pointer. The early pointer specifies a sample value 
which is positioned at least one element before the prompt 
pointer's position, and correspondingly, the late pointer specifies 
a sample value being positioned at least one element after the 
prompt pointer's position. These pointers are used to maintain 

25 the tracking of the received signal by a repeated repositioning of 
the pointers around a correlation maximum value, such that the 
prompt pointer is positioned as close as possible to this value, 
the early pointer is positioned in time somewhat prior to this 
value and the late pointer is positioned in time somewhat after 

30 this value. An important advantage attained by this strategy is 
that one and the same vector may be re-used in order to 
represent different delays. Thus, valuable memory space and 
processing capacity is saved. 

According to still another preferred embodiment of this aspect of 
35 the invention, the tracking involves the following further steps. 
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First, a sequence of incoming level-discrete sample values is 
received. Then, data words are formed from the sample values, 
such that each data word contains a number of elements which 
is equal to the number of elements in each carrier frequency- 
phase candidate vector. After that, a relevant set of carrier 
frequency-phase candidate vectors is calculated for the data 
word, and a pre-generated in-phase representation respective a 
quadrature-phase representation of the vector is acquired for 
each carrier frequency-phase candidate yector in the relevant 
set. Subsequently, each data word is on one hand multiplied 
with a in-phase representation of the carrier frequency-phase 
candidate vector in the relevant set to produce a first inter- 
Tn®d[ate^requency:reduced information word , and on the other 
hand multiplied with a quadrature-phase representation of the 
carrier frequency-phase candidate vector in the relevant set to 
produce a second intermediate-frequency-reduced information 
word. The intermediate-frequency-reduced information words 
thus constitute a respective component of the vectors mentioned 
above that have been derived from the incoming data words. 
This multiplied strategy is advantageous because it enables a 
very efficient usage of the capacity of a general-purpose 
microprocessor. 

According to yet another preferred embodiment of this aspect of 
the invention, XOR- or SIMD-operations are used to perform the 
multiplication between the data word and the in-phase repre- 
sentation of the carrier frequency-phase candidate vector re- 
spective between the data word and the a quadrature-phase 
representation of the carrier frequency-phase candidate vector. 

According to another preferred embodiment of this aspect of the 
invention, the traclcing involves the further steps of: correlating 
the first intermediate-frequency-reduced information word with a 
modified code vector starting at a position indicated by the 
prompt pointer to produce a first prompt-despread symbol string, 
correlating the first intermediate-frequency-reduced information 
word with a modified code vector starting at a position indicated 
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by the early pointer to produce a first early-despread symboi 
string, correlating the first intermedjate-frequency-reduced infor- 
mation word with a modified code vector starting at a position 
indicated by the late pointer to produce a first late-despread 
6 symbol string, correlating the second intermediate-frequency- 
reduced information word with a modified code vector starting at 
a position indicated by the prompt pointer to produce a second 
prompt-despread symbol string, correlating the second inter- 
mediate-frequency-reduced information word with a modified 

10 code vector starting at a position indicated by the early pointer 
to produce a second early-despread symbol string, and corre- 
lating the second intermediate-frequency-reduced information 
word with a modified code vector starting at a position indicated 
by the late pointer to produce a second late-despread symbol 

15 string. Preferably, a resulting data word is then derived for each 
set of the despread symbol strings. This may, for example be 
performed either by performing the relevant adding operation, or 
by looking up a pre-generated value In a table based on the bit 
pattern of the despread symbol strings depending on whether 

20 the resulting data words contain packed or un-packed 
Information. Data of high quality may thus be obtained with a 
minimum of real time calculations, which is advantageous from a 
processing point-of-view. 

According to yet another preferred embodiment of this aspect of 
25 the invention, either XOR- or SIMD-operations are used also to 
perform the multiplication between the intermediate-frequency- 
reduced information words and the modified code vectors, since 
again, this is advantageous from a processing polnt-of-view. 

According to another preferred embodiment of this aspect of the 
30 invention, certain pieces of information are propagated in 
connection with completing the processing of a current data 
word and initiating the processing of a subsequent data word. 
Preferably, this propagated information includes: a pointer which 
indicates a first sample value of the following data word, a group 
35 of parameters that describe the relevant set of carrier 
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frequency-phase candidate vectors, the relevant set of code 
vectors, and the prompt-, early-, and late pointers. Hence, based 
on the propagated Information, the subsequent data word may 
be processed immediately. 

5 According to a further aspect of the invention the object Is 
achieved by a computer program directly loadable into the 
internal memory of a computer, comprising software for perfor- 
ming the above proposed method when said program Is run on a 
computer. 

10 According to another aspect of the invention the object is 
achieved by a computer readable medium, having a program 

recorded- t^^ereonT-iwher^-the -programHs-^o-make- a computer^--:^ 

perform the above proposed method. 

According to another aspect of the invention the object Is 
15 achieved by the signal receiver as initially described, which is 
characterized in that the proposed computer program Is loaded 
Into the memory means, such that the receiver will operate 
according to the above-described method. 

Consequently, the invention offers an excellent Instrument for 
20 receiving and decoding any type of spread spectrum signals In a 
general software radio receiver. For instance, one and the same 
receiver may be used to track navigation satellite signals of 
different formats. These signals may even be received In 
parallel, such that a first signal is received on a first format 
25 while one or more other signals are received on a second or 
third format. 

Moreover, the proposed solution permits one and the same 
receiver to be utilized for completely different purposes, such as 
communication in a cellular system. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

The present Invention Is now to be explained more closely by 
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means of preferred embodiments, which are disclosed as 
examples, and with reference to the attached drawings. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which information may be ex- 
tracted according to the invention, 

Figure 2 shows a spectrum graph over a frequency down- 
converted signal corresponding to the spectrum 
shown In figure 1, 

Figure 3 is an enlarged version of figure 2 in which a 
Doppler shifting of the spectrum is illustrated, 

Figure 4 illustrates how a data signal is modulated onto a 



Figure 5 



signal source specific code sequence according to 
an embodiment of the invention, 

depicts a three-dimensional graph over carrier 
frequency-phase candidate vectors to be used as 
initial correlators according to an embodiment of 
the invention, 

Figure 6 shows an alternative representation the carrier 
frequency-phase candidate vectors in figure 5, 

20 Figures 7-8 show graphs which illustrate how a code index and 
a code phase are defined for a received signal 
according to an embodiment of the invention. 

Figures 9a, b illustrate how a modified code vector is generated 
based on an original code vector according to an 
embodiment of the invention. 



Figure 10 depicts a three-dimensional graph over code rate- 
phase candidate vectors to be used as subsequent 
correlators according to an embodiment of the 

invention, 

30 Figure 11 illustrates how a set of code sequence start 
pointers is assigned to a modified code vector 
according to an embodiment of the invention. 
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illustrates how data words are formed from a 
sequence of incoming level-discrete sample values 
according to an embodiment of the invention, 

illustrates how the data words of figure 12 are 
correlated with various pre-generated vectors 
according to an embodiment of the invention, 

illustrates how multiple pairs of pointers may be 
assigned according to an embodiment of the 
invention, 

shows a proposed signal receiver, and 

illustrates, by means of a flow diagram, the general 
method of jrroces"sTng~¥pre^T"spec^^ 
according to the invention. 

DESCRIPTION OF PREFERRED EMBODIMENTS OF THE 
INVENTION 

In order to realize a truly software based signal processing 
solution in real-time at high sampling rates, such as those 
required in a GNSS, the architecture of the microprocessor 
system must be utilized very efficiently. The present invention 
uses an approach, which involves a form of parallel processing 
based on pre-generated code vector representations. The 
Invention thereby strikes a balance between processing speed 
and memory usage, which compensates for the capacity deficit 
of a general purpose microprocessor in relation to an ASIC. The 
details of this solution will become apparent from the below 
disclosure. 

Figure 1 shows a spectrum graph over an exemplary spread 
spectrum signal from which information may be extracted 
according to the invention. For example, the signal may be a 
GPS C/A-signal in the Li-band with a carrier frequency fnF of 
1,57542 GHz. This means that the signal's spectrum is 



Figure 12 

Figure 13 

Figure 14 

Figure 15 
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symmetrical around the carrier frequency fnp. Moreover, the 
signal has at least 98% of its energy distributed within a 2 IVIHz 
wide frequency band BWhf around the carrier frequency fHp. 
Therefore, it is generally sufficient If the spread spectrum signal 
5 is downconverted to an intermediate frequency fip at around 1 
MHz in order to enable further signal processing at a lower 
frequency. 

Figure 2 shows a spectrum graph over such a frequency down- 
converted signal, whose spectrum is symmetrical around an 
10 intermediate frequency fip, which is considerably lower than the 
carrier frequency fHF, for example fip = 1,25 MHz. Due to move- 
ments of the transmitter (i.e. typically a satellite in orbit) and the 
receiver(i.e; typically a 

may occur in the received signal, which is proportional to the 
15 relative velocity between the transmitter and the receiver. A 
maximum Doppler shift fo is defined which in the GNSS case 
normally is approximately 10 kHz. Thus, the entire spectrum 
may be shifted within an interval 2fD, such that the center 
frequency falls any where between fip+fo (corresponding to a 
20 maximum positive Doppler shift) and fip-fp (corresponding to a 
maximum negative Doppler shift). For illustrative purposes the 
interval 2fD, has been severely exaggerated in the figure 2. The 
figure also indicates a phase variation interval cpo (along an 
independent axis cp) around the intermediate frequency fip. The 
25 phase variation interval (po demonstrates the fact that the phase 
position of the received signal is initially unknown to the 
receiver, and thus constitutes a calibration parameter. The 
effects of this parameter will be further elucidated below with 
reference to figure 5. 

30 An enlarged version of the spectrum in figure 2 is shown in 
figure 3. The intermediate spectrum may, due to possible 
Doppler effects, be shifted such that its center frequency lies in 
the interval 2fD from fiF-min=fiF-fD to fiF-max=fiF+fD- The potential 
maximum shifted spectra are here indicated with a respective 

35 dashed line. The diagram in figure 3 also shows a Nyqulst 
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frequency rs/2 corresponding to a minimal sampling rate rs 
which, when used for sampling the intermediate frequency 
signal, will result in a non-aliased discrete spectrum. Provided 
that fo = 10 kHz, f,F = 1,25 MHz, BWhf = 2 MHz and the inter- 
5 mediate frequency signal is lowpass filtered through an analog 
filter having a 2,5 MHz wide passband, the minimal sampling 
rate rs becomes 5 MHz (i.e. two times the highest frequency 
component after the analog filter). 

Figure 4 illustrates schematically how a data signal D is 
10 modulated onto a signal source specific code sequence CS on 
the transmitter side according to an embodiment of the inven- 
tion. According to a preferred embodiment of the invention, the 

sig rial "sou rce"s^^ 

pseudo random noise. The data signal D here contains a data 
15 symbol sequence [+1, -1, -1, +1, -1], and has a relatively low 
symbol rate, say 50 Hz. The signal source specific code 
sequence CS, however, has a relatively high symbol rate (or 
more correctly chipping rate). For instance, a signal source 
specific code sequence CS in the form of a GPS C/A code may 
20 have a chipping rate of 1,023 MHz and contain 1023 chips per 
period. Each chip, ch, is either +1 or -1. Hence, the C/A code 
repeats itself one per ms. The data signal D is modulated onto 
(or spread by) the signal source specific code sequence CS by 
multiplying each data symbol with the code sequence CS. The 
25 data symbol +1 thereby results in an unaltered code sequence 
CS, whereas the data symbol -1 results in an Inverted code 
sequence CS. 

If a data signal D having a rate of 50 Hz is spread by means of a 
signal source specific code sequence CS having a chipping rate 
30 of 1,023 MHz, this results in 20 entire code sequences CS per 
data symbol. Namely, the period time for one data symbol is 20 
ms, whereas the period time for the code sequence is only 1 ms. 

Figure 5 depicts a three-dimensional graph over different carrier 
frrniirnry phi-r nnrlirhtr irrrtnr . iiiihii h iii iji In iiinrl ni initini < 
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iiii i ii l ii l iii III I iiii lin r i t n in n rn h nr li m mnt ri f tin im ii ii li iiii As 
mentioned above, the intermediate frequency spectrum may be 
Doppler shifted. Therefore, the particular frequency of the 
received signal is initially unknown within a Doppler frequency 
5 interval from fip-min to fiF-max- Moreover, the phase position of the 
received signal is beforehand unknown to the receiver. In order 
to enable an estimation of the actual intermediate carrier fre- 
quency, the Doppler frequency interval fiF-min - fiF-max is divided 
into an integer number of equidistant frequency steps. Of course, 

10 the larger the number of steps is, the more accurately the inter- 
mediate carrier frequency may be estimated. In this example, a 
step size Af of 20 Hz is chosen. In some applications, however, a 

_muc.h_ Ja.rger step sjze Af .JT!ay ._p.rovMe._ a_.s 

accuracy, whereas in other applications (where relatively high 

15 accuracy is needed) a considerably smaller step size Af may be 
required. In any case, a Doppler frequency interval fip-min - fip-max of 
20 kHz and a step size Af of 20 Hz requires 1000 steps, where a 
first step includes the frequencies from f|F-min to f|F.min+Af, a second 
step includes the frequencies from fip-mm+Af to fiF.min+2Af, and so 

20 on up to a last step including the frequencies from fiF-min+999Af to 

flF-max. 

A frequency candidate vector f|F-c is defined for each frequency 
step. In order to estimate an initial phase position of the 
intermediate carrier frequency, an integer number of initial phase 

25 positions <Po, ..., ^7 are determined. The illustrated example 
shows eight such positions <Pc. This gives a phase resolution of 
7i/4. For each frequency candidate vector f|F.c, eight different 
phase positions % ^7 are thus defined which each repre- 
sents a particular phase shift from 0 to 7%/4. A combination of a 

30 certain carrier frequency candidate vector fip-c and a particular 
initial phase position <Pc is referred to as a carrier frequency- 
phase candidate vector. Figure 5 shows all such vectors (i.e. 
eight) for the first step, where fip-c = fiF-min, and the second step, 
where f,F.c = fiF-min+Af. 

35 Additionally, figure 5 includes an axis Lc along which the 
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number of sample values (or length) of the frequency candidate 
vector fiF-c is indicated. The Lc number is preferably a power of 2 
because this generally renders the implementation more simple. 
The number of sample values in each frequency candidate vector 
5 fiF-c is set in consideration of the available memory space, the 
desired frequency resolution and the desired phase resolution. If, 
in the example above, a length of 512 samples is chosen and each 
sample requires one byte, the total (uncompressed) memory 
demand would be 3,9 MB (i.e. 1000 frequencies x 8 phase 
10 positions x 512 B = 4096000 B = 3,9 MB). 

Figure 6 shows an alternative representation the carrier 
frequency-phase candidate vectors Vf<p(fiF.c, ^c) In figure 5. 
HRereT the ve^^^ 
fiF-c=fiF-min (say 1,24 MHz) to fiF.c=fiF-max (say 1,26 MHz), from 
15 the phase shift position 9c=<Po (say 0) to <Pc=^7 (say 7nlA), all 
have a length Lc of Sn samples (say 512) and are represented in 
a cube. Hence, a particular frequency-phase candidate vector 
Vf(p(fiF-c, ^c) is given by a cuboid segment in the cube which is 
parallel with the Lc-axis, and ranges from Lc=Si to Lc=Sn. 

20 After having generated all the carrier frequency-phase candidate 
vectors Vfcp(fiF-c, ^c). the vectors are stored in a digital memory 
with a comparatively short access time, such that they may be 
readily accessed during a later acquisition and/or tracking 
phase. Preferably, the carrier frequency-phase candidate 

25 vectors Vfcp(fiF-c, ^c) are stored according to a data format which 

I is adapted to performing correlation operations between a vooto r m\Mf\MH'f 
that hac h nn n dnrivnri fr n m incoming data words (representing 
one or more received signals) and segmentMDf carrier 
frequency-phase candidate vecto^f(p(f|F-c, ^c)- The vectors may 
ou either include pacl<ed samples or they may contain samples that 
are unpacked to a smallest manageable unit in the micro- 
processor, such as a byte. Further details regarding the data 
format issues will be discussed below. 

Figure 7 shows a graph over an exemplary code sequence CS, 
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which has been received, downconverted and sampled 
according to an embodiment of the invention. After subsequent 
removal of the intermediate frequency component, the code 
sequence CS either attains the signal value +1 or -1, and the 
5 symbols are altered at a chipping rate of for example 1,023 
MHz. At the receiver side, the code sequence CS is sampled 
with a sampling rate rs of for example 5 MHz. The sampling in- 
stances are indicated by means of dots on the code sequence 
signal CS along an axis Sj. As is apparent from the figure there 

10 are approximately four sample values per chip period. Hence, 
the requirements of the Nyquist theorem are clearly fulfilled. 
However, the samples are not capable of determining the 

t!msitimL±n_Majices_wdjJi.j^^ 

code sequence signal CS changes from representing one chip 

15 symbol to representing another one chip symbol. Therefore, a 
so-called code phase must be determined. Moreover, a parti- 
cular sample value at which the code sequence starts must also 
be established. 

Figure 8 shows a graph, which illustrates the sample value at 
20 which the code sequence CS starts (the so-called code index 
CI) and how a code phase, Cph, is defined in relation to this 
sample value. An initial segment of the code sequence signal 
CS in figure 7 is therefore represented in figure 8. 

A first sample instance is here presumed to occur shortly prior 
25 to the actual start of the code sequence CS signal period. 
Hence, a corresponding sample value Si belongs to the end of 
the preceding period. This is illustrated by means of the line 
representing the code sequence signal CS being dashed. A 
subsequent sample instance, however, overlaps the code 
30 sequence CS period. The corresponding sample value S2 is 
defined as the code index CI. The actual beginning of a code 
sequence signal is determined in a preceding acquisition phase 
and typically involves correlation between the received signal 
and a local copy of the code sequence CS. 
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The distance between the start of the code sequence signal CS 
and the code index CI, in turn, is defined as the code phase 
Cph, and thus represents a measure of the translation (or skew) 
between the sample instances Si and the chip transitions. The 
sampling interval between two consecutive sample instances Sj 
is divided into an integer number of possible initial code phase 
positions in order to estimate the code phase Cph. In this 
example, ten 0.1 -steps ranging from 0.0 to 0.9 are defined, 
where a first step 0.0 would indicate that code sequence signal 
CS starts at the code index CI, and a last step 0.9 would 
indicate that code sequence signal CS starts almost at the 
preceding sample value Si. The figure 8 shows a code phase 
Cph of 0. 4. i.e. a situ ation where the code sequence signal CS 
starts approximately in the middle between two consecutive 
sample instances Sj. 

In analogy with the carrier frequency-phase candidate vectors 
Vf(p(fiF.c. ^c), code rate-phase candidate vectors may now be 
defined for each signal source specific code sequence CS and 
combination of code rate CR, code index CI and code phase 
Cph. Nevertheless, according to a preferred embodiment of the 
Invention, a modified code vector is first produced based on 
each original code vector. 

Figures 9a and 9b illustrate how this may be performed. A code 
vector CV in figure 9a contains a number of sample values, say 
5000. Preferably, the number of sample values is selected in 
relation to sampling rate such that a full code sequence period 
is recorded. In the above example, with a sampling rate rs of 5 
MHz, this means that the 5000 samples cover exactly one full 
code sequence CS period (which has a duration of 1 ms and 
contains 1023 chips, since the chipping rate is 1,023 MHz). 

A modified code vector CVm is now generated on basis of the 
original code vector CV by copying a particular number of 
elements Ee (e.g. two) from the end of the original code vector 
CV to the beginning of the modified code vector CV^. 
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Correspondingly, a particular number of elements Eb (e.g. two) 
from the beginning of the original code vector CV are copied to 
the end of the code vector CVm, i.e. after the elements Ee. 
Figure 9b shows the resulting modified code vector CVm. Since, 
5 in this example, each chip symbol is sampled at approximately 
four instances, two sample values are equivalent to about one 
half chip period. Thereby, the above copying of the elements Eg 
and Eb allows a one-half-chip shifting of the code vector CV 
during a correlation operation between the code vector CV and 
10 the received data. This will be discussed further below with 
reference to figure 1 1 . 

Figure 10 depicts a three-dimensional graph over different code 
rate-phase cancR 

a signal source specific code sequence CS(i). Analogous to the 
15 intermediate carrier frequency, the nominal code rate CRc for 
the code rate candidate vector CRc-c may also vary due to 
Doppler effects from a lowest value CRmm to a highest value 
CRmax- However, since the code rate (or chipping rate) is com- 
paratively high, the rate variation here becomes rather small. 

20 For example, a Doppler shift of 10 kHz at the intermediate 
carrier frequency is equivalent to a shift in a nominal code rate 
CRc at 1,023 MHz of merely 6,5 Hz. Namely, the Doppler shift in 
respect of the code rate is generally equal to the ratio between 
the carrier frequency and the code rate (1575,42 MHz / 1,023 

25 MHz = 1540; a Doppler shift of 10 kHz at the intermediate 
carrier frequency thus corresponds to 10000/1540 « 6,5 Hz). 
The code rate candidate vector CRc-c may attain a value any- 
where in a Doppler rate interval between CRmm = 1022993,5 Hz 
and CRmax = 102306,5 Hz. The Doppler rate interval CRmin - 

30 CRmax is divided into an integer number of equidistant code rate 
steps ACR, say 13. In this example, each interval CRmin, 
CRmin+ACR, CRmax thus spans 1 Hz. According to a preferred 
embodiment of the invention, the code rate step size is adaptive, 
e.g. determined on basis of the sampling frequency and the 

35 desired resolution. 
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In figure 10, the three-dimensional graph (or cube) contains a set 
of code rate-phase candidate vectors CVm(CRc-c. Cph), where the 
code rate ranges from CRc-c=CRmin to CRc.c=CRmax. the code 
phase ranges from Cph=0.0 to Cph=0.9, and each vector has a 
5 length Lcvm of Sm samples (for example 5004). Hence, a particular 
code rate-phase candidate vector CVm(CRc.c, Cph) is given by a 
cuboid segment CVm(CRmin+2ACR, 0.4) which is parallel with the 
Lcvm -axis and ranges from Lcvm=Si to Lcvm =Sm- 

After having generated all the code rate-phase candidate 
10 vectors CVm(CRc.c, Cph), these vectors are stored in a digital 
memory with a comparatively short access time, such that they 
may be readily accessed during the acquisition and/or tracking 
'phase. Prefer^bry, fhe^cWe rate- phase "can did a fe~~Vecf6rs' 
CVm(CRc-c, Cph) are stored according to a data format which is 
15 adapted to correlation operations to be performed between a 
vector that has been derived from incoming data words and a 
segment of a the code rate-phase candidate vectors CVm(CRc-c, 
Cph). For example, the code rate-phase candidate vectors 
CVm(CRc-c, Cph) may include packed samples or they may 
20 contain samples that are unpacked to a smallest manageable 
unit in the microprocessor, such as a byte. Further details 
regarding these data format issues will be discussed below with 
reference to figure 13. 

During tracking of a spread spectrum signal, i.e. a continued 
25 reception of the signal and demodulation of the data contained 
therein, it is necessary to update relevant tracking parameters 
for the signal. Normally, the tracking phase is preceded by a so- 
called acquisition phase during which a set of preliminary 
parameters are established that are required for initiating the 
30 signal decoding. A successful acquisition phase thus results in 
that at least one signal source specific code sequence is 
identified and a data signal transmitted by means of this 
sequence is possible to demodulate. The tracking characteris- 
tics that may be associated with a signal source specific code 
35 sequence include: a modified code vector CVm (which defines a 
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carrier frequency candidate vector fiF.c and an initial pliase 
position 9c), a code pliase position, Cpli, and a code index, CI. 
In order to maintain an adequate timing of the correlating 
operations in the receiver, a set of pointers are required which 
5 indicate the start position of the code sequence in relation to the 
modified code vector and are updated repeatedly, preferably 
between each correlation. 

Figure 11 illustrates one such set of pointers that are assigned 
to a modified code vector CV^. A prompt pointer Pp indicates a 

10 current estimate of the code sequence start position. This 
pointer Pp is Initially set equal to the code index CI. Adequate 
Pp-polnter positions for any subsequent segments of the 

^ received' data are obtained by retrieving appropriate modified 
code vectors CVm from the digital memory where such pre- 

15 generated are stored. Between each segment, the code rate 
candidate vector CRc-c as well as the initial code phase Cph is 
updated. Additionally, at least one pair of early- and late 
pointers Pe and Pl respectively is assigned on each side of the 
prompt pointer Pp, where the early pointer Pe specifies a sample 

20 value being positioned at least one element before the prompt 
pointer's Pp position, and the late pointer Pl specifies a sample 
value being positioned at least one element after the prompt 
pointer's Pp position. Each correlation is then performed with the 
presumption that the code sequence starts at the Pp-position, at 

25 the PE-position as well as at the PL-position. If the correlation 
over the Pp-position results in a higher correlation value than 
over any of the other positions, and at the same time the Pe- and 
PL-pointers are balanced (i.e. the correlation over these 
positions result in equal values), this is interpreted as an optimal 

30 setting of the tracking parameters. Should, however, this not be 
true the sampled data must be repositioned in relation to the Pp-, 
Pe- and PL-pointers, such that the correlation peak again 
coincides with the Pp-pointer's position. An un-calibrated pointer 
positioning may thus be detected via a difference between the 

35 correlation value over the PE-ppinter position and the correlation 
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value over the PL-pointer position, wliile the correlation over the 
Pp-posltion still results in the highest correlation value. In this 
case, the potentially highest correlation value (corresponding to 
the optimal Pp-positioning) lies somewhere between the current 
5 Pp-position and the current PE-position or between the current 
Pp-posltion and the current Pt-position depending on which of 
the Pe- and Pt-pointer positions that results in the highest 
correlation value. A pointer updating normally takes place 
between each segment. A second or third pair of Pe- and Pl- 

10 pointers outside the first pair of pointers further enhances the 
possibilities of accomplishing an accurate adjustment of the 
tracking. Depending on the value of the received symbol, an 
optimal setting Pp-pointer may, in fact, also be equivalent to the 
correlation over the Pp-position resulting in a lower correlation 

15 value than over any of the other pointer positions. Namely, this 
Is true if a negative valued data bit is received. However 
otherwise, the same principles apply. 

In any case, the Pp-, Pe- and Pt-pointers may only be positioned 
within certain intervals in the modified code vector CVm. The Pp- 

20 pointer may be set anywhere within a first interval 111, the Pe- 
pointer within a second interval 112 and the Pt-pointer within a 
third interval 113. All the pointers are always set together, such 
that their relative distances remain unchanged. A limit value PL-min 
for the PL-pointer typically coincides with the first element of the 

25 modified code vector CVm. Namely, outside this range a 
meaningful correlation cannot be performed. As a consequence of 
this and the fact that the pointers are set jointly Pp-, Pe- and Pl, a 
corresponding limit value PE-min for the PE-pointer is also defined. 

Figure 12 illustrates how data words d(1), d(2), d(N) are 
30 formed from a sequence 1210 of incoming level-discrete sampJe 
values according to an embodiment of the invention. The data 

words d(1) d(N) are formed such that each data word d(k) 

contains a number of elements which is equal to the number of 
elements Sn in each carrier frequency-phase candidate vector. For 
35 instance, the length of these vectors may be 512 sample values 
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(as in the example described above witli reference to the figures 
5 and 6). Given that each sequence 1210 contains 5000 sample 
values, the number N of data words becomes 10 (5000/51 2»9, 76 
means that the last quarter of the tenth data word d(10) is empty). 
5 According to the invention, the elements in a current data word 
d(k) are processed jointly (i.e. in parallel). Then, the processing of 
a subsequent data word d(k+1) is initiated, and any processing 
parameters obtained during the processing of the current data 
word d(lc) that are required when processing the subsequent data 

10 word d(k+1) are propagated as input data to the processing of the 
latter. These processing parameters preferably include: a pointer 
Pd which indicates a first sample value of the subsequent data 
word d(k+1), a group of parameters which describe the relevant 
set of carrier frequency-phase candidate vectors Vf^pCfip-c, ^c) (i e. 

15 a carrier frequency candidate vector fiF-c representing an in-phase 
version and a carrier frequency candidate vector fip.c representing 
a quadrature-phase version), the relevant set of code vectors 
CVm, and prompt-, early-, and late pointers Pp, Pe, and Pl 
respectively. 

20 Figure 13 illustrates how the data words of figure 12 are corre- 
lated with various pre-generated vectors according to an embo- 
diment of the invention. First, a relevant set of carrier frequency- 
phase candidate vectors Vf,p(f|F.c, ^c) is calculated for the data 
word d(k). This set is normally obtained In a preceding 

25 acquisition phase (or, if the data word d(k) is not the first word 
to be processed, during the processing of a previous word). 

Thereafter, for each carrier frequency-phase candidate vector 
Vf(p(fiF-c, ^c) in the relevant set, a pre-generated in-phase repre- 
sentation f|Fi and a quadrature-phase representation fipQ of the 

30 vector {Vf^{fif.c, ^c)) respectively is acquired. Both these repre- 
sentations are found in the digital memory where the pre-gene- 
rated carrier frequency-phase candidate vectors \/f^{f\f.c, ^c) are 
stored because they are merely one quarter period delays of 
one and the same vector, and the memory holds initial phase 

35 position representations for an entire period, see figure 6. 
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Then, in order to eliminate tlie influence of the carrier frequency 
component, the data word d(k) Is multipiied with the in-phase 
representation fipi of the carrier frequency-phase candidate 
vector Vf<p(fiF-c. *Pc) in the relevant set. As a result thereof, a first 
5 intermediate-frequency-reduced information word S|F-i(k) Is pro- 
duced. Correspondingly, the data word d(k) is also multiplied 
with the quadrature-phase representation fipQ of the carrier 
frequency-phase candidate vector Vf<p(f|F-c, ^c) in the relevant 
set, and a second intermediate-frequency-reduced information 
10 word S|F-Q(k) is produced. 

Preferably, the multiplication between the data word d(k) and 
the in-phase representation f|Fi respective the quadrature-phase 
representation fiFQ of the carrier frequency-phase candidate 
vector Vf<p(fiF-c, ^c) are performed by means of XOR- or SIMD- 
15 operations (i.e. 1-bit and multiple-bits multiplications respec- 
tively). This is namely possible if the data word d(k) and the 
vector Vf(p(fiF-c, ^c) have compatible data formats. 

Subsequently, the first intermediate-frequency-reduced informa- 
tion word SiF-i(k) is multiplied with a modified code vector CVm-p, 

20 which has been retrieved from the digital memory where these 
pre-generated vectors are stored, and starts at a position 
indicated by the prompt pointer Pp. A first prompt-despread 
symbol string A|p(k) is produced as a result of this operation. 
The first intermediate-frequency-reduced information word 

25 SiF-i(k) is also multiplied with a modified code vector CVm-E(l<), 
which starts at a position indicated by the early pointer Pe. A 
first eariy-despread symbol string AiE(k) is produced as a result 
of this operation. The first intermediate-frequency-reduced 
information word SiF.i(k) is likewise multiplied with a modified 

30 code vector CVm.L(k), which starts at a position indicated by the 
late pointer Pl, and a first late-despread symbol string AiL(k) is 
produced. Moreover, the second intermediate-frequency- 
reduced information word SiF-Q(k) is multiplied with a modified 
code vector CVm.p(k), which starts at a position indicated by the 

35 prompt pointer Pp, and a second prompt-despread symbol string 
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AQp(k) is produced. Correspondingly, the second Intermediate- 
frequency-reduced Information word S|F.Q(k) Is multiplied with a 
modified code vector CVn,.E(k), which starts at a position 
indicated by the early pointer Pe, and a second early-despread 
5 symbol string AQE(k) Is produced. Likewise, the second 
Intermediate-frequency-reduced information word S|F-Q(k) is 
multiplied with a modified code vector CVm-Uk), which starts at a 
position indicated by the late pointer Pl, and a second late- 
despread symbol string AQL(k) is produced. 

10 After that, a resulting data word DR.iE(k), DR.ip(k), DR.iL(k), 
DR.QE(k), DR.Qp(k) and DR.QL(k) may be derived from the each of 
the despread symbol strings A|p(k), A|E(k), A|L(k), AQp(k), AQE(k) 
and AQL(k) respectively by adding the elements in the respective 
string together. If the strings represent un-packed data, the 

15 processor may simply perform the relevant adding operation to 
obtain the resulting data word DR(k). If however, the strings 
represent packed data, the resulting data word DR(k) is, 
according to a preferred embodiment of the invention, derived by 
looking up a pre-generated value in a table 1310 based on the 

20 bit pattern given by the respective despread symbol string 
A,p(k), A,E(k), A,L(k), AQp(k), AQE(k) and AoUk). 

According to a preferred embodiment of the invention, in addition 
to the individual data words DR.iE(k), DR.ip(k), DR.iL(k), DR.QE(k), 
DR-Qp(k) and DR.QL(k), corresponding accumulated data words are 
25 produced, such that after having generated a data word DR(k), the 
sum of all data words Dr(1) to DR(k) is also obtained. Hence, 
when the data word d(N) has been processed, the resulting sum 
Dr(1) to Dr(N) is aiso at hand. 

Moreover, the resulting data words represent payioad information. 
30 For instance, based on the set of resulting data words DR.|E(k), 
DR.ip(k), DR.iL(k), DR.QE(k), DR.Qp(k) and DR.QL(k) a demodulated 
piece of information of a transmitted data symbol sequence D as 
shown In figure 4 may be derived. In a steady state operation, 
basically only one of the resulting data words, say DR.ip(k), 
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actually carries the information. 

According to a preferred embodiment of the invention, the multi- 
plications involved in the multiplications betvi^een the interme- 
diate-frequency-reduced information v/ords SiF-i(k) and S|F-Q(k) 
5 respectively and the modified code vectors CVm-p(k), CVni-E(k); 
CVm-L(k) are all performed by means of XOR- or SIMD-ope- 
rations. This is possible if the intermediate frequency infor- 
mation words SiF-i(k); S|F-Q(k) and the modified code vectors 
CVm-E(k), CVm.p(k); CVm-L(k) have compatible data formats. 

10 It was mentioned above with reference to figure 11 that more 
than one pair of early- and late pointers may be assigned to the 
prompt pointer in order to accomplish an improved tracking. 
Figure 14 illustrates how two pairs of such pointers Pei; Pli and 
Pe2; Pl2 respectively are assigned around a prompt pointer Pp 

15 according to an embodiment of the invention. The horizontal 
axis here represents a code shift CS and the vertical axis shows 
a normalised correlation factor. An optimal positioning of the 
prompt pointer Pp is thus equivalent to a normalised correlation 
value of 1. The more the sampled data is shifted in either 

20 direction from this position, the lower the correlation value 
becomes down to substantially zero for a shift of one chip or 
more. The algorithm is perfectly balanced when (as illustrated in 
the figure) both the pointers Pei and Pli of a first pair produce a 
first shifted correlation result, both the pointers Pe2 and Pl2 of a 

25 second pair produce a second shifted correlation result, the first 
shifted correlation result is larger than the second shifted corre- 
lation result and both the first and the second shifted correlation 
results are smaller than the correlation result obtained at the 
prompt position Pp. Said correlation results are preferably 

30 generated according to a process corresponding to what has 
been described above with reference to figure 13. Based on the 
correlation results for the pointer positions Pp, Pei, Pli, Pe2 and 
Pl2 a (if necessary new) set of modified code vectors CVm.E(k), 
CVm.p(k); CVm-L(k) is selected, which is deemed to be optimal. 
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Figure 15 shows a signal receiver 1500 for receiving navigation 
data signals transmitted in a navigation satellite system accor- 
ding to an embodiment of the invention. The receiver 1500 
includes radio front end unit 1510, an interface unit 1520 and a 
6 digital processor unit 1530. 

The radio front end unit 1510 is adapted to receive a continuous 
radio signal Shf. and In response thereto produce a corre- 
sponding electrical signal Sif Which has a comparatively high 
frequency. The interface unit 1520 is adapted to receive the 

10 electrical signal Sif, and in response thereto, produce a 
sequence of sample values that represents the same information 
as the electrical signal S|f and is divided into data words d(k). 
The digital processor unit 1630 is adapted to receive the data 
words d(k), and in response thereto, demodulate a data signal. 

15 The digital processor unit 1530, in turn, includes a memory 
means 1535, which is loaded with a computer program that is 
capable of controlling the steps of the proposed procedure when 
the program is run in the processor unit 1530. 

In order to sum up, the general method of processing spread 
20 spectrum signals according to the invention will now be 
described with reference to a flow diagram in figure 16. 

A preparation step 1600 pre-generates code vectors which each 
represents a signal source specific code sequence that is 
intended to be received (or at least possible to receive) and 
25 demodulated in the receiver. The step 1600 is performed before 
the signal reception is Initiated. 

A step 1610 then receives a continuous signal of a compara- 
tively high frequency. A following step 1620 samples the 
continuous signal at a basic sampling rate, whereby a resulting 
30 sequence of time discrete signal samples is produced. Each 
sample is also quantised (either with a relatively low-resolution, 
such as with 1 bit per sample value, or with a relatively high 
resolution depending on the application and the number of data 
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bits per samples delivered from the radio front end unit 1510), 
such that a corresponding level-discrete sample value is 
obtained. Subsequently, a step 1630 forms data words from the 
sample values, where each data word includes one or more 
5 consecutive sample values. After that, a step 1640 carries out 
correlation operations between the information In the data words 
and at least one of the pre-generated representations of a signal 
source specific code sequence. The correlation step Involves 
correlating at least each vector in a sub-group of the code 

10 vectors with at least one vector that has. been derived from a 
data word. A following step 1650 produces data as a result of 
the correlation performed in step 1640. Then, a step 1660 

.investigates ._whetheL 

whether there are any more data words to process. If more 

15 sampled data to process is found, the procedure returns to the 
step 1640. Otherwise, the procedure loops back to the step 
1610 again for continued reception of the incoming signal. 
Naturally, such signal is preferably also received during 
execution of the steps 1620 to 1660. The sequential procedure 

20 illustrated in figure 16 is merely applicable to a specific received 
signal segment. Preferably, all the steps are actually 1610 to 
1660 are actually performed in parallel. 

The process steps, as well as any sub-sequence of steps, 
described with reference to the figure 16 above may be 

25 controlled by means of a programmed computer apparatus, such 
as a microprocessor located In a GNSS-receiver. Moreover, 
although the embodiments of the Invention described above with 
reference to the drawings comprise computer apparatus and 
processes performed in computer apparatus, the invention thus 

30 also extends to computer programs, particularly computer 
programs on or in a carrier, adapted for putting the invention 
into practice. The program may be in the form of source code, 
object code, a code intermediate source and object code such 
as in partially compiled form, or in any other form suitable for 

35 use in the implementation of the process according to the 
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invention. Tlie carrier may be any entity or device capable of 
carrying the program. For example, tlie carrier may comprise a 
storage medium, such as a ROM (Read Only Memory), for 
example a CD (Compact Disc) or a semiconductor ROM, or a 
5 magnetic recording medium, for example a floppy disc or hard 
disc. Further, the carrier may be a transmissible carrier such as 
an electrical or optical signal which may be conveyed via 
electrical or optical cable or by radio or by other means. When 
the program is embodied in a signal which may be conveyed 
10 directly by a cable or other device or means, the carrier may be 
constituted by such cable or device or means. Alternatively, the 
carrier may be an integrated circuit in which the program is 
em bedded , the.integrated_.drcuit_b pe rform i ng , 

or for use in the performance of, the relevant processes. 

15 The term "comprises/comprising" when used in this specification 
is taken to specify the presence of stated features, integers, 
steps or components. However, the term does not preclude the 
presence or addition of one or more additional features, 
integers, steps or components or groups thereof. 

20 The invention is not restricted to the described embodiments in 
the figures, but may be varied freely withjn the scope of the 
claims. 
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Claims 

1. A method of processing spread spectrum signals com- 
prising: 

receiving a continuous signal (S|f) of a comparatively high 
5 frequency, 

sampling the continuous signal (Sif) at a basic sampling rate 
(rs), whereby a resulting sequence of time discrete signal 
samples (S[Si]) is produced, 

quantising each signal sample into a corresponding level- 
10 discrete sample value, 

forming of a plurality of data words (d(1), .... d(N)) which 
each includes one or more consecutive sample values (Si, Sn), 
and 

correlation between information in the data words (d(k)) and 
15 at least one representation (CS(i)) of a signal source specific 
code sequence (CS), characterized by 

the method comprising a preparation for the correlation step, 
wherein, before receiving the continuous signal (Sif), a multitude 
of code vectors (CV; CVm) are pre-generated, each code vector 

20 representing a particular code sequence (CS(i)) of the at least 
one signal source specific code sequence (CS), and 

the correlation step involving multiplying at least each vector 
in a sub-group (CVm-E, CVm-p; CVm-i) of the code vectors (CV; 
CVn,) with at least one vector (SiF.|(k); SiF.Q(k)) derived from the 

25 data word (d(k)). 

2. A method according to claim 1, characterized by each 
code vector (CV; CV^) representing a particular signal source 
specific code sequence (CS(i)) being sampled at the basic sam- 
pling rate (rs) and quantised with the quantising process being 

30 used to produce the level-discrete sample values (si, .... Sn). 

3. A method according to any one of the preceding claims, 
characterized by the at least one signal source specific code 
sequence (CS) being pseudo random noise. 
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4. A method according to any one of the preceding claims, 
characterized by the receiving step involving down conversion 
of an incoming high-frequency signal (Shf) to an Intermediate 
frequency signal (Sif), the high-frequency signal (Shf) having a 
5 spectrum which is symmetric around a first frequency (fHp) and 
the intermediate frequency signal (Sif) having a spectrum which 
Is symmetric around a second frequency (fiF) being considerably 
lower than the first frequency (fMp). 



5. A method according to claim 4, characterized by the steps 
10 of: determining a maximum frequency variation (fp) of the 

second frequency (fip) due to Doppler-effects, 

defining a Doppler frequency interval (fiF.min - fip-max) around 

the second frequency (fip), the Doppler frequency interval (f|F.min - 

fiF-max) having a lowest frequency limit (fiF-min) equal to the 
15 difference between the second frequency (fip) and the maximum 

frequency variation (fp), and a highest frequency limit (fiF-max) equal 

to the sum of the second frequency (fip) and the maximum 

frequency variation (fo), 

dividing the Doppler frequency interval (f|F-min - fiF-max) into an 
20 integer number of equidistant (Af) frequency steps (fiF-min, fiF-min+Af, 

■ fiF-max), and 

defining a frequency candidate vector (fiF-c) for each fre- 
quency step (fiF.min, fiF-min+Af, fiF-max). 



6. A method according to claim 5, characterized by the steps 
25 of: determining an integer number of initial phase positions 
(^0, ^7) for the intormodiato frcquGnay o i gnal (Srrt , and 

defining a carrier frequency-phase candidate vectoi 
(Vfq,(fiF-c, ^c)) for each combination of carrier frequency 
candidate vector (fiF-c) and Initial phase position (9o, <P- 



30 7. A method according to claim 6, characterized by the steps 
of: determining the number of elements (Sn) in each carrier 
frequency-phase candidate vector (Vf,p(fiF.c,^c)), and 
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storing the carrier frequency-phase candidate vector 
(Vf(p(fiF-c, ^c)) according to a data format being adapted to per- 
I forming correiatipn operations between 4 h9 at l ejbl uiiy vuutor mHiplfcd 
\ 4 fiiT..|(i i ), : ii.u (IO) ilLi i ueU fiui i i t he data word and a 

5 segment of a carrier frequency-phase candidate vector (Vf<p(fiF.c, 
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8. A method according to claim 7, characterized by the 
adapting of the data format involving adding at least one ele- 
ment to each segment of the carrier frequency-phase candidate 
vector (Vf<p(f|F-c, ^c)) such that #ve segment attains a-number of 
elements which is equal to the number of elements in ■ each th e 

(d(k)), thus enabling a segment of the carrier frequency-phase 
candidate vector (Vfq,(f|F-c, ^c)) and one of the at least one 
vector (SiF.i(k); Sif-q(I<)) to be processed jointly by at least one 
of a SIMD-operation and an XOR-operation. 



9. A method according to any one of the preceding claims, 
characterized by the steps of: 

determining a maximum variation a code rate due to 
20 Doppler-effects, 

defining a Doppler rate interval (CRc-min - CRc-max) around 
a center code rate (CRc), the Doppler frequency interval having 
a lowest code rate limit (CRc-min) equal to the difference 
between the center code rate (CRc) and the maximum code rate 
25 variation (CRd), and a highest frequency limit (CRc-max) equal to 
the sum of the center code rate (CRc) and the maximum code 
rate variation (CRd), 

dividing the Doppler rate interval (CRc-min - CRc-max) into 
an integer number of equidistant (ACR) code rate steps (CRc-min, 
30 CRc-min+ACR, .... CRc-max) and 

defining a code rate candidate (CRc-c) for each code rate 

step (CRc-min. CRc-min+ACR, .... CRc-max)- 
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10. A method according to claim 9, characterized by deter- 
mining an Integer number of possible Initial code phase positions 
(0.0, 0.9) for each code rate candidate (CRc-c). 

11. A method according to claim 10, characterized by 
5 defining, for each signal source specific code sequence (CS(i)), 

a set of combinations between code rate candidate (CRc-c) and 
code phase position (0.0, 0.9), each combination repre- 
senting a code rate-phase candidate vector. 

12. A method according to claim 11, characterized by 
10 generating, for each signal source specific code sequence 

(CS(i)), a set of code vectors (CV) by- 
sampling each code rate-phase candidate vector with the 
basic sampling rate (rs) whereby a corresponding code vector 
(CV) is produced. 

15 13. A method according to any one of the preceding claims, 
characterized by generating a modified code vector (CVm) on 
basis of each code vector (CV) by: 

copying a particular number of elements (Eg) from the end 
of an original code vector (CV) to the beginning of the modified 
20 code vector (CV^), and 

copying the particular number of elements (Eb) from the 
beginning of the original code vector (CV) to the end of the code 
vector (CVm). 

14. A method according to claim 13, characterized by storing, 
25 for each signal source specific code sequence (CS(i)), a set of 

modified code vectors ({CVm(CRc-c, Cph)}), where 

each modified code vector (CV^) contains a number of 

elements (Si, Sm) representing a sampled version of at least 

one full code sequence (CS), 
30 a particular modified code vector (CVm) is defined for each 

combination of code rate candidate (CRc-c) and code phase 
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position (Cph). 

15. A method according to claim 14, characterized by 
adapting tiie data format of the modified code vectors (CVm) with 
respect to the data format of the at least one vector (S|F-i(k); 

5 SiF.Q(k)) derived from the data word {6{k)), such that a modified 
code vector (CVm) and one of the at least one vector (S|F.|(k); 
Sif-q(I^)) may be processed jointly by at least one of a SIMD- 
operation and an XOR-operation. 

16. A method according to any one of the claims 14 or 15, 
10 characterized by involving an initial acquisition phase and a 

subsequent tracking phase/ wherein duHng 
a set of preliminary parameters are established which are re- 
quired for initiating a decoding of signals received during the 
tracking phase, a successful acquisition phase resulting in at 
15 least one signal source specific code sequence (CS) being identi- 
fied and a transmitted signal being rendered possible to track, 
each of the at least one signal source specific code sequence 
(CS) being associated with tracking characteristics in the form of: 

a modified code vector (CVm), 
20 a carrier frequency candidate vector (fiF-c), 

an initial phase position (9c), 

a code phase position (Cph), and 

a code index (CI) denoting a starting sample value for the 
modified code vector (CVm). 

25 17. A method according to claim 16, characterized by the 
tracking involving: 

calculating, based the tracking characteristics, a prompt 
pointer (Pp) for each modified code vector (CVm), the prompt 
pointer (Pp) indicating a code sequence start position, and an 
30 initial prompt pointer (Pp) being equal to the code index (CI), and 
assigning, around each prompt pointer (Pp), at least one 
pair of early- and late pointers (Pe, Pl; Pel Pli Pe2. Pl2), where 
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the early pointer (Pe) specifies a sample value being positioned 
at least one element before the prompt pointer's (Pp) position, 
and the late pointer (Pl) specifies a sample value being 
positioned at least one element after the prompt pointer's (Pp) 
5 position. 

18. A method according to claim 17, characterized by the 
tracking involving: 

receiving a sequence of incoming level-discrete sample 
values (1210), 

10 forming data words (d(1), d(N)) of the sample values 

(1210) such that each data word (d(k)) contains a number of 
elements which is equal to the number of elements (Sn) in each 
carrier frequency-phase candidate vector (yf^ihf-c, ^c)). 

calculating a relevant set of carrier frequency-phase candi- 

15 date vectors (Vf,p(f|F-c, ^c)) for the data word (d(k)), and 

acquiring, for each carrier frequency-phase candidate 
vector (Vf<p(f|F-c, ^c)) in the relevant set, a pre-generated in- 
phase representation (fipi) and a quadrature-phase represen- 
tation (fiFo) of the vector (Vfp(f|F.c, ^c)) respectively. 

20 19. A method according to claim 18, characterized by the 
tracking involving: 

multiplying each data word (d(k)) with the in-phase repre- 
sentation (fiFi) of the carrier frequency-phase candidate vector 
(Vf(p(fiF-c, ^c)) in the relevant set to produce a first Intermediate- 

25 frequency-reduced information word (S|F-i(k)), 

multiplying each data word (d(k)) with the quadrature- 
phase representation (fiFo) of the carrier frequency-phase candi- 
date vector (Vfq,(fiF-c, ^c)) in the relevant set to produce a second 
intermediate-frequency-reduced information word (Sif Q(k)). 

30 20. A method according to claim 19, characterized by the 
tracking involving: 

multiplying the first intermediate-frequency-reduced infer- 
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mation word (S|F-i(k)) with a modified code vector (CVm-p(i<)) 
starting at a position indicated by tlie prompt pointer (Pp) to 
produce a first prompt-despread symbol string (A|p(l<)), 

multiplying the first intermediate-frequency-reduced infor- 
5 mation word (SiF.i(k)) with a modified code vector (CVm-E(k)) 
starting at a position indicated by an early pointer (Pe) to 
produce a first early-despread symbol string (A|E(k)), 

multiplying the first intermediate-frequency-reduced infor- 
mation word (SiF.i(k)) with a modified code vector (CVm-L(k)) 
10 starting at a position indicated by a late pointer (Pl) to produce 
a first iate-despread symbol string (A|L(k)), 

multiplying the second intermediate-frequency-reduced in- 

fprmatipn word (SiF..Q(k)) with a. modified code vector (CVm.p(k)) 

starting at a position indicated by the prompt pointer (Pp) to 
15 produce a second prompt-despread symbol string (AQp(k)), 

multiplying the second intermediate-frequency-reduced in- 
formation word (S|F-Q(k)) with a modified code vector (CVm-E(k)) 
starting at a position indicated by the early pointer (Pe) to 
produce a second early-despread symbol string (AQE(k)), and 
20 multiplying the second intermediate-frequency-reduced in- 

formation word (S|F-Q(k)) with a modified code vector (CVm.L(k)) 
starting at a position indicated by the late pointer (Pl) to 
produce a second Iate-despread symbol string (AQL(k)). 

21. A method according to claim 20, characterized by the 
25 tracking involving deriving, for each despread symbol string, a 

resulting data word (DR.ip(k), DR.iE(k), DR.iL(k), DR.Qp(k), DR.QE(k); 

DR.QL(k)). 

22. A method according to claim 21, characterized by deriving 
the resulting data words (DR.ip(k), DR.iE(k), DR.iL(k), DR.Qp(k), Dr. 

30 QE(k); DR.QL(k)) by looking up a respective pre-generated value in 
a table (1310). 
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23. A method according to any one of the claims 19 - 22, 
characterized by performing the multiplication between the data 
word (d(k)) and the in-phase representation (fipi) of the carrier 
frequency-phase candidate vector (Vfq,(fiF-c, ^c)) respective 

5 between the data word (d(k)) and the a quadrature-phase repre- 
sentation (fiFo) of the carrier frequency-phase candidate vector 
(Vfcp(fiF-c, ^c)) by means of at least one of a SIMD-operation and 
an XOR-operation 

24. A method according to any one of the claims 19 - 23, 
10 characterized by performing the multiplication between the inter- 
mediate-frequency-reduced information words (SiF.i(k); SiF.Q(k)) 
and the modified code vectors (CVm.p, CVm-n; CVm.L) by means of 
at least one of a SIMD-operation and an XOR-operation. 

25. A method according to any one of the claims 19 - 24, 
15 characterized by propagating, in connection with completing 

the processing of a current data word (d(k)) and initiating the 

processing of a subsequent data word (d(k+1)): 

a pointer (pd) indicating a first sample value of the 

subsequent data word (d(k+1)), 
20 a group of parameters describing the relevant set of carrier 

frequency-phase candidate vectors (Vf,,(fiF-c, ^c)), 
the relevant set of code vectors (CVm), and 
prompt-, early-, and late pointers (Pp, Pe, Pl)- 

26. A computer program directly loadable into the internal 
25 memory of a computer, comprising software for controlling the 

steps of any of the claims 1 - 25 when said program is run on 
the computer. 

27. A computer readable medium, having a program recorded 
thereon, where the program is to make a computer control the 

30 steps of any of the claims 1 - 25. 
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28. A signal receiver (1500) for receiving navigation data 
signals transmitted in a navigation satellite system comprising: 

a radio front end unit (1510) adapted to receive a conti- 
nuous radio signal (Shf) and in response thereto produce a 
5 corresponding electrical signal (S|f) comparatively high 
frequency, 

an interface unit (1520) adapted to receive the electrical 
signal (S|f) and In response thereto produce a sequence of 
sample values being divided Into data words (d(k)), and 
10 a digital processor unit (1530) adapted to receive the data 

words (d(k)) and In response thereto demodulate a data signal, 
and Including a memory means (1535), characterized in that 
a computer program according to clairn 26 Is loaded Into the 
memory means (1535). 

15 
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Abstract 

The present invention relates to processing of spread spectrum 
signals, where a continuous signal of a comparatively high 
frequency is received (1610). This signal is sampled at a basic 
5 sampling rate whereby a resulting sequence of time discrete 
signal samples is produced, which are in turn quantised into a 
corresponding level-discrete sample value (1620). A plurality of 
data words are formed, which each includes one or more 
consecutive sample values (1630). Information obtained from 

10 these data words is correlated (1640) with at least one 
representation of a signal source specific code sequence, which 
has been pre-generated in the form of a code vector (1600). The 
correlatibn step^ specifically cbn^lMing at'Teas^^^ 

vector in a sub-group of the code vectors with at least one 

15 vector that has been derived from the data word. Thereby 
resulting data is produced (1650). 
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